Achievable Qubit Rates for Quantum Information Wires 
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Suppose Alice and Bob have access to two separated regions, respectively, of a system of electrons 
moving in the presence of a regular one-dimensional lattice of binding atoms. We consider the 
problem of communicating as much quantum information, as measured by the qubit rate, through 
this quantum information wire as possible. We describe a protocol whereby Alice and Bob can 
achieve a qubit rate for these systems which is proportional to N~'S qubits per unit time, where 
N is the number of lattice sites. Our protocol also functions equally in the presence of interactions 
modelled via the t — J and Hubbard models. 
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One of the key requirements for a functioning quantum 
information processor is the ability to transport quantum 
information from one location to another (see, e.g., [ij 
for further details.) Hence, finding physical systems us- 
able as error-tolerant quantum communication channels 
has become a priority. Recently a class of physical sys- 
tems realising such channels, namely strongly interacting 
quantum spin systems, has been investigated [3-[l)@"16|. 
Under fairly general conditions it has been shown that 
such quantum information wires can communicate quan- 
tum information with arbitrarily high fidelity. 

Generally speaking, the proposals developed so far ex- 
plore hilbert space rather inefficiently: typically, once a 
signal has been placed in the wire, Alice must wait until 
Bob has successfully removed the signal before another 
signal can be sent. In other words, the qubit rate Q — 
the number of qubits which can be successfully communi- 
cated per unit time — for these protocols scales roughly 
inversely with the length of the system, i.e., Q ~ N~^, 
where N is the number of lattice sites. So far there has 
only been one proposal 0^ exploring qubit rates be- 
yond this: there a protocol for an engineered chain is 
described which achieves, after taking energy normalisa- 
tion into account, a qubit rate Q ^ . 

In this Letter we aim to improve this situation in 
two ways: firstly, we study naturally occurring one- 
dimensional systems of particles moving in the presence 
of a regular lattice of binding atoms, and secondly de- 
scribe a protocol achieving a qubit rate Q N~3. Our 
results do not require that the interactions between the 
particles are engineered, however, we do assume that the 
sender and receiver each have good control of at most 
one lattice site each. 

We initially study free fermions moving in the pres- 
ence of a regular one-dimensional lattice of N atoms. 
The number iV 3> 1 is considered to be a large pa- 
rameter which is divisible by 4 for convenience. We 
denote by {fljl^Li the operators annihilating a fermion 
at site j; the operators {aj}jLi obey the canonical 



anticommutation relations {ajjaj.} = 6jk- Thus our 
hilbert space is fock space, and has an orthonormal basis 
given by the occupation number basis: |ni, 77.2, riAr) = 

(al)"'(4)"' •■• (a]v)""l^)' where nj e {0,1} are the 
number of electrons located at site j, and we identify 
the vacuum state \Q) = |0, 0, . . . , 0). We initially model 
the lattice geometry using the tight-binding model [l7| . 
describing free electrons hopping on a regular lattice of 
N sites: 
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a^flj+i + h.c. 



(1) 



where we identify site (A^ + 1) with site 1, i.e., we assume 
periodic boundary conditions. The periodic geometry is 
a theoretical device and we later argue that our results 
extend straightforwardly to the more realistic chain set- 
ting. We also argue later that our protocol functions 
equally well in the presence of interactions modelled by 
the t — J model and the Hubbard model. 

We define the length of the system to be 1. 
Therefore, the single-particle quantum state |j) = 
|0, 0, . . . , Ij , . . . , 0) means that there is a particle sitting 
with probability one at position Xj — j /N in the lattice. 
We refer to the subspace T-Ls spanned by |j) as the single- 
particle subspace. When we refer to \j) we say that the 
particle is at site j corresponding to physical position Xj . 

As the number of lattice sites increases, and if the 
quantum state of the system is the discretisation of a 
sufficiently smooth wavefunction, the dynamics will be 
equivalent to that of a free particle on a circle. (See the 
supplementary material for further details). The con- 
sequence of this is that gaussian-modulated wavepack- 
ets with wavenumber k propagate at some group velocity 
v{k) with essentially no change in shape. 

In our protocol Alice wants to communicate a sequence 
of (possibly unknown) qubit states lipa), en = 1,2,..., M, 
to Bob. Alice and Bob are each allowed access to a region 
Ra - {1,2,..., vN"^, {Rb = {7V/2, N/2 + 1,..., N/2 + 
vN^ — 1}, respectively) , of vN^ sites, where is a con- 
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stant to be chosen later, at opposite ends of the ring. 
(We later argue that Alice and Bob only need access 
to one site within their regions, respectively.) Initially 
the system is prepared in the vacuum state jfi). Alice 
and Bob each have access to M ancilla qubits labelled 
Aa, (respectively, Be): the qubits Aa are initialised in a 
product of the M qubit states |V-'a) and Ba arc initialised 
in a product of some convenient fiducial state |0). Thus 
the initial state of the system before the protocol begins 

is \A)M---\MAMmO)B^---\0)BM- 

To accomplish the communication task, Alice performs 
some encoding operation U\ coupling her region Ra and 
the qubit Ai containing the first message IV-"!)- The sys- 
tem is now in the state t/i|V'i)|f^)- (From now on, when 
discussing the state of the system we suppress mention of 
the qubits which haven't yet interacted with the lattice.) 
The system evolves as \^{t)) = e-'"'^Ui\\pi)\Vt) . After 
a time T, where T depends only on H and UN's, Bob 
performs a decoding operation Vi on his addressable sites 
Rb in order to decode or refocus the communicated state 
into the ancilla qubit B\. The state of the system at the 
moment Bob has applied Vi is Vie-^"'^Ui\i]ji)\n)\Q). The 
protocol is deemed to succeed when the average fidelity 
of the decoded state p\ = trg^(Vie~'^'^?7i|V'i)(V'i| 8' 
O \0){0\ule'"'^V^), where tr^. denotes the partial 
trace over all the system+ancillas except for Bj, given by 

Tr^^JdQ,{i>,\pM, (2) 

where \tpi) is the input state averaged over the Bloch 
sphere, is above some prespecified threshold value 1 — 
e. In the case where Alice wants to send M signals we 
require that each of the average fidelities J^a, of the M 
messages is above 1 — e. We shall later obtain a lower 
bound for the average fidelity of the protocol we're about 
to describe. 

The operations Ua are designed to simply swap an un- 
known qubit state Itpa) from register a into her accessibil- 
ity region Ra; they are realised by applying the unitary 
operator 

where cr± = ^{a^ T i^^l), {(^a^cry,a^} are the 
Pauli sigma operators on qubit A^, and g = 

iSjei?^ ^ ^"^ ^'^ ''^3^ with (T, I, and k to be cho- 
sen later. The number 7 is chosen so as to ensure that g 
obeys the canonical anticommutation relation {g,g^} = 

1, i-e.,{ff,ff^} = 7^ Eje-RA ^ = 1- Note the useful 
identity 

f^a = I - CT+CT-fif^t _ a-a+g^g + a+g + a'g^. (4) 

The operator g^ has been chosen so as to create a fermion 

in a discretisation of a gaussian envelope centred on Z = 
Nxq in Alice's region with wavenumber k. 



To understand how our encoding operation works we 
write IV'i) = ci|0)-|-di|l). The transformation 
i7i|Vi)|0) isgivenby C/i|V'i)|n) = ci|0)|O)+rfi|0)5t|fi) = 
|0)(cil -I- dig^)\ri), it simply leaves the |0) state alone 
and flips the state |1) and creates a single fermion in a 
gaussian- modulated wavefunction in Ra- 

We can now understand the dynamics of the encoded 
state |*(0)) = |0)(cil + dig^)\n): we find that 

|*(t)) = |0)(ciI + di5Hi))|O), (5) 

where g{t) = e~'^^* ge^^^ . This is, in turn, given by the 
single-particle dynamics g{t) = X]j j'=i .9j(0)['5**^]jj'aj', 

where 5(0) = g, gj{Q) = 76" ^+2^*'=^ j & Ra, and 
(j'j(O) = otherwise, and A is the N x N matrix whose 
matrix elements are given by = Sjj'+i + 6jj'-i, 

l<j,f<N. 

The operation Ui corresponds to depositing a fermion 
into the lattice in a discretisation of a gaussian modulated 
wavefunction. After some time t the particle will have 
propagated out of the region Ra- Thus, Alice's region 
will be indistinguishable from the vacuum jQ). In this 
case the action of the next encoding operation U2 will be 
proceed just as for Ui- However, the application of U2 
may not proceed due to some residual amplitude for the 
particle to remain in region Ra and consequently to be 
swapped back to the ancilla A2 . Such a situation counts 
as an error. We now argue that the magnitude of this 
error is captured by the quantity 

\\{g{0),g\t)}\\ = \{gmgm, (6) 

where \g{t)) = g{t)\Cl) = J2j 9j{t)\j) is a single-particle 
state. Let's temporarily assume that this quantity is zero 
for all t- In this case there is no amplitude for a particle 
to be swapped back into the ancillas and the state of the 
system after M messages are created is given by 

|*m) = Umb-'^'Um-i ■ ■ ■ C/2e-'^*C/i|Vi) • • • IVm)]^^) 
= |0)(cmI + W(0)) • • • (cil + dig^iM - l)t))\n), 

(7) 

In this case we can intepret the state |^'m) as that of 
M free fermions moving independently without any in- 
terference effects whatsoever. 

After a certain time T has elapsed the first signal will 
have reached Bob's decoding region. Bob then decodes 
the message into his first ancilla by applying the opera- 
tion 

_ ^-i^{a+a-hh^+a-a-^h^h)^i^{a+h+a-h^)^ ^g-j 

where acts on qubit and h = j'^^j^n^ 93^)0'h 
where 7' is chosen to ensure the anticommutation rela- 
tion {/i, /?,^} = 1. If wc temporarily assume that .g(T) is 
given by g{T) — j'^j^j^^ gj{T)aj we see that Vi is the 
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optimal way to decode the signal as it completely swaps 
the fermion out of the lattice into the first ancilla. Evi- 
dently in the case where \g{T)) does not have this form 
then the decoding operation is not completely successful. 
Again a key role is played by the quantity 



broadening L{t) due to dispersion coming from the cubic 
correction we use the third- order broadening factor [20| 



\{g{T),hm\^MT)\h)\, 



(9) 



when it is equal to 1 the decoding is completely success- 
ful. This quantity, when not unity, is interpreted as the 
discretisation error eo- 

Let's now consider the realistic situation where 
II{5(0)j5^(0}II 7^ 0; it is clear that the system state is 
no longer of the form |^'m) given by Eq. ([7]). In general 
we have that 



■ f/2e^'"*f/i|Vi) • • • \i^M)\n) - |*M> + ITm), 

(10) 

where now Iv^m) is of the form Eq. ([7]) (but no longer as- 
sumed normalised) and ITm) is an orthogonal state which 
we interpret as an error. Using a combination of the tri- 
angle inequality, the leibniz property of the commutator, 
the anticommutation relations for g and g^ , and the fact 
that < 1, we obtain the bound 

eE = \\\Tm)\\ < 3 ^ {M - j)\{gmgUt))\. (11) 
j=i 

where e^; is the encoding error. Thus the errors oc- 
curring during the running of our protocol can be com- 
pletely understood in terms of the overlaps {g{0)\g{t)) 
and {g{t)\h). Both of these quantities pertain to the 
single-particle sector; we have reduced our problem to 
understanding the dynamics of a single particle propa- 
gating through the lattice. This problem is now very 
well understood in the literature, and a comprehensive 
study may be found, for example, in Q. We summarise 
the pertinent results here (see also the supplementary 
material). Firstly, the gaussian modulated wavefunction 
\g{0)) with wavenumber k propagates through the sys- 
tem at the group velocity v{k) — '^^}^'' = — ^ sin {^k) 
through the system. Choosing fcg = N/A maximises this 
group velocity and we obtain v{ko) = —Att/N (recall that 
Alice is located at physical position and Bob's at po- 
sition 1/2). Thus it takes a time t — N/{8Tr) for the 
wavepacket to travel from Ra to Rb- Secondly, as long 
as the wavepacket is broad enough it retains its gaussian 
shape (up to some small errors) throughout the propa- 
gation from Alice's region to Bob's region. It turns out 
that an initial width ~ A^3 is sufficient to ensure that the 
wavepacket doesn't disperse too quickly. This is because 
a wavepacket of physical width L(0) corresponds, in mo- 
mentum space, to a wavepacket of width L~^(0) centred 
on wavenumber k — N/A. Around this wavenumber the 
dispersion relation uj{k) = 2cos(^fc) is approximately 
linear, with a small cubic correction. To estimate the 
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(12) 



We solve this equation for t ~ N/ (Stt) and find that 
a choice of L(0) ^ Na is sufficient to ensure that the 
wavepacket spreads by at most a constant amount. Ob- 
viously there are errors coming from the truncation of 
the Alice's wavepacket to the region Ra, and from the 
discretisation. One can check that all these errors can be 
made exponentially small by redefining the constants a 
and \Ra\ (these calculations are presented in the supple- 
mentary material). 

We find \{g{0)\g{t))\ g-'^^t^'N 3^ ^here k is a con- 
stant that only depends (logarithmically) on the desired 
errors ep and eE, where ep is the propagation error. 
Thus, if we want to reduce the total encoding error aris- 
ing from |||rjv/)|| to eE it is sufficient to wait a time 
t ^ vNSj where the constant v can be chosen to scale 
as v ^ — log(e) . The decoding error experienced by Bob 
is proportional to \{h\g{T))\, where T ^ iV is the time 
taken for a wavepacket to enter Bob's region Rb- As has 
been argued previously Q, this quantity is simply re- 
lated to the weight of the wavefunction in Bob's region. 
By choosing \Rb\ to scale as TVs sites this quantity can 
be increased towards 1 exponentially quickly. 

The average fidelity J^a for each of the decoded qubits 
Ba is thus straightforwardly bounded from below by 
J- a > 1 — eE — ep — ep). As we've argued, these three 
sources of error are, in turn, reduced exponentially fast 
in c by redefining a, \Ra\, and by a constant factor 
c. 

It is inevitable that during the running of the protocol 
errors will build up in the system. Thus the achievable 
rate will also be implicitly determined by the maximum 
running time Tmax desired because the errors depend on 
M , the total number of messages to be sent: the error 
estimate above only applies if M < iVs. In the case 
where M = XN's we need to add up the errors linearly 
in A > 1 in the obvious way. Once the accumulated error 
goes above some prespecified threshold it is necessary to 
cool the system down again to the ground state. If we 
assume this takes some constant time then we learn that 
the achievable rate is not affected. 

To conclude the description of our protocol we have 
to argue that Alice and Bob can carry out the encod- 
ing and decoding steps with access to only a single site. 
It turns out that a method to achieve this has already 
been described in Q in the spin-system setting. One can 
check that it applies with the obvious modifications to 
the situation we have here. 

Our protocol applies to the tight-binding model of elec- 
trons propagating around a ring. However it also extends 
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without modification to the chain geometry. Ah that 
needs to be checked is that the relevant results concerning 
the free propagation of waves extend to the line-segment 
geometry. This is a straightforward exercise and is left 
to the reader. 

There are two natural extensions of our model. The 
first is to the setting where fermion-fermion interactions 
are included and the second is to the setting where sig- 
nals are encoded in different initial wavepackets. In the 
first generalisation we model interactions via the (spin- 
less) t — J model: Ht-j — —tJ2^=i'^]^j+i + ^.c. -I- 
J'^J^i^ rijUj+i. (Our analysis also applies to the case 
with spin.) In this case we may apply the previous analy- 
sis by noticing that ||e-^-«"'-'"|5'A/)-e-*-f^^|^'M) || < \s\ei. 

This follows from triangle inequality and not- 
ing that the particle-particle interaction term 
Hj = "^J^i^ njUj^i, when applied to I^'a/), satis- 
lies \\Hi\<SJM)\\^ei. 

The magnitude of e/ may be reduced exponentially by 
linearly increasing t, the time between signals. This is 
because the state of the system is comprised of well- 
separated single particles, and hence particle-particle in- 
teractions are negligible. A similar argument applies to 
the Hubbard model i?HM = -tY^j.ai^l,cr'^j+i,o- + h.c.) -I- 
[/ ^^-^ nj.j.nj_-|-, where Cj^a annihilates an electron with 
spin a at site j. In this case, our protocol proceeds as be- 
fore after identifying Uj with Cj^^ and ignoring the other 
spin degree of freedom. Because of the onsite interac- 
tions it is impossible to exploit the second spin degree of 
freedom to run two instances of the protocol in parallel 
without some nontrivial modification. 

The second generalisation is for Alice to create addi- 
tional signal fermions in wavepackets ik{x) orthogonal 
to the gaussian wavepacket. Indeed, the hermite poly- 
nomials naturally suggest themselves here. The nearly 
linear dispersion relation suggests that these orthogo- 
nal wavepackets also propagate without change of shape, 
and hence provide additional communication channels 
between Alice and Bob. However, a hermite polynomial 
is very close to a translation of a gaussian. This implies 
that these wavepackets propagate at a different group 
velocity. If one attempts to find N-^^^ such orthogonal 
wavepackets one runs into the problem that the approx- 
imations involved in the derivation of the group velocity 
and dispersion relation begin to break down, especially at 
the higher frequencies. This generalisation may provide 
a constant qubit rate but is unlikely to work for many 
realistic models (except for the tight-binding model) due 
to the fermion-fermion interactions. 

We have described a protocol whereby Alice and Bob, 
having access to only small regions of a system of elec- 
trons moving in the presence of a lattice of binding atoms, 
can communicate quantum information with arbitrarily 
high fidelity at a rate of N^'s qubits per unit time. This 
result improves considerably upon the extant protocols 



for quantum information wires in two important ways. 
Firstly, the rate we achieve here is much greater than 
that encountered in the quantum spin-system case, and 
improves on the rate achieved by [6| for a specially engi- 
neered chain. Secondly, our protocol applies to naturally 
occurring systems modelled by the tight-binding, t — J, 
and Hubbard models. There is considerable scope for fur- 
ther work following our contribution: the extension of our 
analysis to the bosonic case, e.g., to the Bose-Hubbard 
model modelling cold atoms in optical lattices, certainly 
merits investigation. Also, a more detailed analysis of the 
case where different encoding wavefunctions may yield 
better — possibly even constant — qubit rates. 
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Supplementary material 
Bounding the error terms 

We bound the error term \Tm) as follows: using the identity Eq. (j4]) we have that 

UMe''"'---e-'"'Ui\i^i)---\i^Mm) = |*M) + |rA/), (13) 

where 

IFm+i) = |Ai) + IA2) + IA3) + IA4), (14) 

and 

|Ai) =[/M|VM+i)e-*^*|rAf-i), 

IA2) =dM|l)(I-ff(0)<?t(0))e-«*|*M-i), 

|A3)--CM|0)gt(0)g(0)e-'^*|*M-i>, 

IA4) =CM|l)5(0)e-^^*|vI'M-i). 

The norms of each of these four terms can be bounded as follows. Firstly we have by unitary invariance of the norm 
that |l|Ai)|| < |||rM-i)||- The second may be bounded by noticing that 

|||A2)|| < \dM\\M0)9\0), (cM-iI + rfA/-i.9^(l)) • • • (cil + d.g^m - l)t)m)\\ 

< E M^ IIK^M-iI + dM-igH^)) ■ ■ ■ (c.+il + d,+^g\{M - J - mMO)g\0),g\{M - j)t)]x 

(c,_il + d,_igt((M - J + l)t)) ■ ■ ■ (cil + digt((Af - l)t))|r!)|| 

< E Um.gHiM - jm\ \\{cM-il + dM-wHl)) ■ ■ ■ (c,+il + d,+,g\{M ~ J - l)i))5^(0)x 

(cj_il + d,-ig^((M - J + l)t)) ■ ■ ■ (cil + dig^iM - l)t))|17)|| 

AI-l 

< E ll{.9(0),5^((M - j)t)}\\ ||gt(0)|| ||(cM-iI - dM-igH^)) ■ ■ ■ (c,+il - d,+^g\{M - j - l)t))x 
3=1 

(cj_il + d,-ig^((A/ - J + l)t)) ■ ■ ■ (cil + dig\iM - l)t))\n)\\ 

< E \\{9{0),gH3t)}l 

where we've used the triangle inequality, the leibniz property of the commutator, the fact that \\g\\ < 1, and that 
IKca/I + '^A/3^(l)) ■ • • (C2I + dig'' {Mt))\n)\\ < 1 for all Cj and dj. The third term is bounded using exactly the same 
argument, and we obtain 

Af-l 

|||A3)||< E Kff(o)l3(J*))l- (17) 

The final term is bounded using anticommutation relations for g(0) and g{t): 

M-l 

g{0)e-^"'\^M-i) = E d,{gmgm - j)t)){cM-il- dM-ig\l)) ■ ■ ■ {c,+,l ^ d,+,g\{M - j - l)t))x 

(c,_il + d,-ig\{M - j + l)t)) ■ ■ ■ (cil + dig\{M - l)t))\^l). (18) 
Taking the norm gives us, via the triangle inequality, 

Af-l 

|||A4)||< E K3(0)I3(J*))I- (19) 



(16) 
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Thus we learn that 

M-1 M-1 

\\\rM)\\ < ll|rA^-i)|| + 3 5] MOMm < 3 {m ~ j)\{g{0Mjt))\. (20) 



Dynamics in the single-particle sector 

In this appendix we review some resuhs concerning the propagation of discrete gaussian wavepackets for scalar 
particles propagating in a regular one-dimensional lattice. 

Throughout this appendix we work in the single-particle subspace T-Ls- This is the subspace of hilbert space with 
basis \ j) = aj|r2), j = 1, 2, . . . , TV. Since the tight-binding model Eq. ([1]) preserves particle number then if the system 
begins in "Hs it remains there for all time. The first thing to note is that the matrix-elements of Eq. ^ in T-Ls are 
given by 

[A],fe = {j\H\k). (21) 

where A is the N x N matrix whose matrix elements are given by [AJ^j' = ^j.j'+i + Sjj'-i, 1 < < N, where we 
identify the first site with the {N + l)th site: 1 ee -|- 1. 

Since A is circulant it may be diagonalised via a discrete fourier transform; we obtain 

N 

A = Y,uj{k)\W{k)}{W{k)\, (22) 
fe=i 



where 



uj{k) = 2 cos [^k] (23) 



1 ^ 

|iy(fc))^^^M^'^-|j), (24) 



is the dispersion relation with corresponding eigenvectors 

N 

and where /i is the TVth root of unity, ^ = e^^*. 

As noted in the body of the paper, the operator A may be understood as the discretised kinetic energy operator. 
Hence, the dynamics generated by A are physically equivalent to that of a freely propagating particle. These are now 
familiar facts and, with the appropriate qualifications, may be made mathematically rigourous (see, e.g., any book 
on the numerical analysis of PDE, say). For simplicity, rather than summarise the relevant (large) mathematically 
rigourous literature, we content ourselves here with a discussion at the level of physical rigour. The reader may be 
assured that all statements here can be brought to the level of mathematical rigour with a suitable amount of technical 
effort. 

The first step is extend, via the natural map, every state |0) = Y^^=i 4'{^)\W {k)) of our discrete system to a 
hmction : ^ C in L2(S^): 



N 

|,^) ^(/.(x) ==^0(/c)efe(x), (25) 

k=l 

where ek{x) — e^'^^''^ and x G [0, 1). Note that (j) is normahsed with respect to the standard L2 inner product. In this 
way we see that 

m = cj^ix,), (26) 
where Xj = j/N. Recall that the momentum operator p = —i^ acts on </> as a multiplication operator 

N 

p(f)^2n'Yk4>{k)ek{x). (27) 

k=l 
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Now the image A of the operator A acts as the multipUcation operator 



N 



A<p = J2^ik)mek{x), (28) 
fe=i 

where ui{k) = 2 cos {^k). We now focus on momenta k near ko — ^ (assuming 4| A?") to maximise the group velocity: 
suppose that 

k = ko+l, (29) 
with I smaU in comparison to ko. We then expand the dispersion relation w(A;) 

.{k) = .(ko) + l-iko) + ^^(fco) + ^^(fco) + • • • 



(30) 

0- ^1 + 0+ + ■■■ . 



7v' + ^+3! 7V3 



(j){x) = (j){l)eko+iix) dl = Ckoix) (j){l)ei{x)dl = eko{x)(l)o{x), (31) 



Now suppose that 



where A is a cutoff. Clearly, for A o{N), we have that the action of A on is well approximated by the action of 
the wave operator 

^~->+3!l^^" 

on (po{x), i.e., 

Acf>{x) = Bk, (x) [-^P+ + (33) 

In particular, with the choice A = nNi , we have that the action of A on 0o(a;) is equivalent to 

^M^) ~ + M^) + 0{N-i). (34) 

Thus, as long as the discrete fourier representation of \<p) obeys 

N feo+A 

\cl,) = ^^{k)\W{k))= m\W{k)), (35) 

fe=l fc=fco— A 

we can approximate the dynamics of A via the dynamics of the wave operator —jfP + -g^^P^- Evidently this 
approximation will only hold for t ~ o{Ni). However, this is more than enough for our situation as t ~ A'' is the time 
it takes for the wavepacket to traverse the distance between Alice and Bob. 

Now we investigate the smallest width in real space a wavepacket obeying the constraint A = k;A''3 can have. To 
do this we introduce the dirac comb 

r?iv(a;) = ^^<5(x-jyiV). (36) 

Suppose that there is a L2(S^) function 0(a;) such that (/)(x) = for a; ^ (0,1) and = 4'{xj) (indeed, such a 

function always exists, and further, cj) may be chosen to be C°°). In this case we have that 
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for k = 1, 2, . . . , iV, where 

/■oo 

nfi^m) = / e'-'''-f{x)dx, (38) 



is the fourier transform operator. In this way we see that 

N 



^J2i^*f}m\wik)). (39) 



Now we recall that a gaussian 



k=l 



<?(x) = (40) 



has fourier transform 



g{k) = 7r3V2^e2"'^^«e"2'^ . (41) 



The objective is to choose (t){x) — g{x). Unfortunately such a choice doesn't satisfy (f>{x) = for x ^ (0, 1). (And, 
indeed, such a wavefunction is nontrivial outside Alice's region.) However, we can choose a so that this constraint is 
satisfied up to an exponentially small error. 

Let's now calculate the a required to approximately satisfy the constraint (j){x) — for x ^ Ra and the constraint 
4>o{k) = for fc ^ A] required by Eq. psp . The second constraint is satisfied up to an exponentially small error 
e~'^ with the choice 

= (42) 

This corresponds, in real space, to a gaussian 17(2;) with characteristic width 

m = (43) 

This corresponds, in lattice units, to a region Ra of width \Ra\ = NL{Q) = sites. 

Now, under the dynamics generated by A a gaussian wavefunction does not remain a gaussian. Instead, it becomes 
an Airy function [20j. To estimate the width of the evolved gaussian \g{0)) ~ e~**^|(7(0)) (which is found from 
g{x,t) = e~'^*^g{x,0)) we need to use the third-order broadening factor (20l | 



m 
m 



1 + 7: 



1 /w"'(fco)t^ ^' 



2 V\/2L3(0) 



(44) 



We solve this equation for t = iV (the time it takes to traverse the lattice): we find that with our choice of that L(0) 
such a wavepacket spreads by at most a constant amount. 

The third-order broadening factor may derived by estimating the overlap |(<7(0)|(7(i)) | for arbitrary t. Using the 
fact that 

(g(0)|g(i)) = (5(0),e-"^5(0)) (45) 
we reduce our problem, via Parseval's relation, to the estimation of the Fourier- Airy integral: 

2<jV^ r e-4-'-''=%^*^-tt^*'=' dk. (46) 

J —00 

Choosing t — ^xiN^ (the time it takes for a wavepacket to propagate from Ahee's region) we have that 

pOO 2 8 
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Notice that, because of the gaussian weighting, the integrand is nontrivial only for k < kNu , so that changing variables 

to fc = kN^I we have 



We now expand in the smallncss of the third exponent: 

^ e-2-'%2"'«'^i (^1 - ^{27rfN-ixiKH^ + ---^dk. (49) 

The dominant term is evidently the first term in the expansion (this is only valid as long as Xi is 0(1), i.e., on 
timescales longer than the lattice scale; for very short times Xi ~ -/V~^ the expansion breaks down): 

e (50) 

Higher order terms come from differentiation of e 2<= with respect to k, and the remainder R{k) is easily estimated 
to scale as 

= const. X e ^ , (51) 

where the constant depends only on c and k. Recall that xq < xi < , where Xg corresponds to the centre of Alice's 
region and Xi = Ns /2 corresponds to a separated physical position. By increasing xi we can reduce this overlap to 
zero doubly exponentially fast. Hence we have our estimate 

l(ff(0)|ff(i)>l~const. xe-^^. (52) 

This tells us that to reduce the total error 

\\\rM)\\<\{9mgm (53) 

to less than e we need only wait a time t xN^ between signals. 



